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Abstract 



We investigate the self-organizing nature of relativistic quantum field theory 
in terms of canonical transformation and duality presenting simple but ex- 
plicit examples of (4> 4 )i + i and (0 6 )i + i theories. Our purpose is fulfilled by 
applying the oscillator representation (OR) method which allows us to convert 
the original strong interaction theory into a weekly interacting quasiparticle 
theory that is equivalent to the original theory. We discuss advantages of the 
OR method and compare the results with what was already obtained by the 
method of Gaussian effective potential (GEP) and the Hartree approximation 
(HA). While we confirm that the GEP results are identical to the Hartree re- 
sults for the ground state energy, we found that the OR method gives the 
quasiparticle mass different from the GEP and HA results. In our exam- 
ples, the self-organizing nature is revealed by the vacuum energy density that 
gets lowered when the quasiparticles are formed. In the (0 6 ) n~i theory, we 
found two physically meaningful duality-related quasiparticle solutions which 
have different symmetry properties under the transition of quasiparticle field 
However, these two quasiparticle solutions yield the identical effec- 
tive potential in the strong coupling limit of the original theory. 
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I. INTRODUCTION 



With the current advances of Relativistic Heavy Ion Collision(RHIC) physics, there is 
a growing interest in discussing the self-organizing nature of relativistic quantum field the- 
ory (RQFT). The phase transition and the spontaneous symmetry breaking anticipated to 
be observed in the RHIC facilities are the paramount examples of the physical phenomena 
due to the self-organization of quantum fields. Since these novel phenomena cannot be eas- 
ily predicted in the ordinary perturbation series, they form highly nontrivial examples that 
ought to be analyzed by various nonperturbative methods available for the RQFTs. 

Although the problem of quantization of interacting fields is not yet completely solved, 
it appears that the duality existing in the RQFTs is the key to handle the strong inter- 
action theories. According to the duality, the original particle theory of strong couplings 
may be equivalent to the quasi-particle theory of weak couplings. As the interactions get 
stronger, the associate quantum fluctuations get larger and consequently the formation of 
a nontrivial vacuum often accompanied by the condensation of fields may be energetically 
more favorable than the adherence to an original trivial vacuum without any condensation. 
The particles moving in the ambient nontrivial vacuum condensates may be described by the 
quasi-particles which carry different masses compare to their original masses defined in the 
trivial vacuum. Then, the couplings among the quasi-particles may become weak as the cou- 
plings among the original particles grow strong. When such duality works self-consistently 
within the RQFT, one may be able to solve the corresponding strong interaction problem 
utilizing the weak coupling developed for the quasi-particles. For the phenomena involv- 
ing only weak interactions, the perturbation theory may provide systematic predictions and 
even some physical insight how the interaction can affect the properties of the system. Thus, 
the duality allows us to convert the original nonperturbative problem into the quasiparticle 
perturbative problem. However, the key issue here is to check if the RQFTs indeed have 
the self-organizing characteristics that lead to the duality within themselves. To investigate 
such nontrivial characteristics, one cannot just rely on the perturbation series from the be- 
ginning but need to develop a nonperurbative method which may be useful to analyze the 
self-organizing nature of the RQFTs. 

In this work, we discuss such a nonperturbative method that appears to be much simpler 
than other rather well-known methods such as the Hartree Approximation (HA) and the 
Gaussian Effective Potential (GEP) method. Following the nomenclature in the literature 
we call this method as the Oscillation Representation (OR) method, although in our view all 
of these nonperturbative methods (OR,GEP,HA) stem from the same principle of quantum 
effective action approach [1]. The OR method was explicitly formulated by Efimov [2] and 
was, in fact, used earlier by Chang [3] and Magruder [4], though in some indirect way. The 
method is at length described in a monograph [5]. The basic idea of the OR method is to 
redefine the mass of interacting field relative to the free one and simultaneously introduce a 
shift of the field quantization point leading to a nonzero value of its vacuum condensate. This 
effect is realized in the nature represented by a spontaneous symmetry breaking mechanism 
and thus one can expect that it may be a general feature of quantum field systems with 
interactions. However, the OR method has not been utilized as extensively as the methods 
of GEP [6-10] and HA [11,12]. Even in the 1+1 dimensional scalar field theory such as the 
one we present in this work, not much of the OR results distinguished from the GEP or HA 
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results have been known or discussed to the best of our knowledge. Although the OR results 
for the (0 4 )i + i theory were presented in Ref. [5], the distinction between the OR and GEP 
(or HA) results for this simple model was not pointed out at length. Thus, for the present 
work, we analyze the details of the distinction in this most simple case first and extend the 
application to the (0 6 )i+i theory: 

C = ^0^0 - im 2 2 - - H\ (1) 

where the results from the variational method of the GEP have also been studied [9]. The 
GEP is defined as 

VgOo) = min n < <f> , Q\H\Q, O >, (2) 

where 7i is the Hamiltonian density corresponding to Eq.(l). Here, Q is the mass parameter 
and \ fl, 0o > is a normalized Gaussian wave functional centered on = O , i.e. 

< O , O >= 1, < 0o, fi|0|Q, 00 >= 00- (3) 

As we will discuss the details in the next two sections (Sections II and III), the GEP method 
can in principle give the two OR equations that we solve simultaneously in this work. One 
of them is known as the f2 equation in the GEP,z.e. d < O , f2|7Y|f2, O > /dVL = at 
Q = Cl, and the other OR equation is equivalent to c?Vg(0o)/^0o — 0. However, in the GEP 
method, these two equations are not solved simultaneously in practice. It turns out that 
the quasiparticle mass defined by d 2 Vc(4>o) /d<fro m the GEP is not same with which is the 
simultaneous solution of the quasiparticle mass from the two OR equations. In this paper, 
we therefore pay more attention to the distinct aspects of the OR results compared to the 
GEP and HA results and discuss the property of the OR results which was not mentioned 
in the more great detail before. 

Formally, the OR is settled by requiring that the Hamiltonian operator should be written 
in terms of creation and annihilation operators of an oscillator basis with an appropriate 
frequency and in the correct form defined by 

(1) all field operators in the total Hamiltonian H = Hq + Hj are written in the normal 
ordered product, 

(2) the Hamiltonian H is quadratic over the field operators, 

(3) the interaction Hamiltonian Hj contains field operators in powers more than two. 
Using this method with the cannonical transformation of quantum fields, one may allow 
to build a standard perturbation theory for the quasiparticles in the region of the coupling 
constants where the original theory is nonperturbative. Thus, we utilize the OR method to 
transform the original theory which is intrinsically nonperturbative to the equivalent theory 
of the quasiparticle representation that can be solved rather easily by the perturbation 
method. 

As we will discuss, the exact OR results include several vacuum solutions yielding dif- 
ferent mass values of the quantum field for fixed values of coupling constants. Among those 
multiple solutions, however, it is certainly not difficult to select the non-trivial solutions 
that are physically meaningful because of the duality consideration. In this work, we will 
thus focus only those duality-related solutions which are not hindered from the physical 
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interpretation. Nevertheless, we note that the OR method generates also the solutions that 
cannot be interpreted as straightforward as the ones we focus in this paper. In particular, 
we find a class of solution that doesn't exhibit the anticipated duality property and another 
class of solution that we call spurious in the sense that it yields the pure imaginary value for 
the vacuum condensation. We do not yet know all the physical meaning of these classes of 
solutions, but the generation of multiple solutions are expected because the OR equations 
are the nonlinear algebraic equations relating the quasiparticle mass and the field conden- 
sation with the parameter of coupling constant as we will present in details for the rest of 
the paper. 

In Section II, we present more details of the OR method using the 4 theory which can 
be easily deduced from the Lagrangian density of c/> 6 theory by taking h = 0. We com- 
pare the OR results with the results obtained by the other available methods such as the 
GEP and the HA as well as the available lattice result. In Section III, we apply the OR 
method to the (0 6 )i+i theory and discuss the effects from the presence of 6 interaction to 
the duality-related vacuum solutions. We analyze the quasiparticle mass and the energy 
density in the entire domain of coupling parameter space and present the classical potential 
for the nontrivial solutions. We discuss two physically meaningful duality-related quasipar- 
ticle solutions with different symmetry properties and show that they however generate the 
identical effective potential in the strong coupling limit of the original theory. Summary and 
conclusions follow in Section IV. In the Appendix A, we briefly summarize some details of 
the HA in the (0 4 )i + i theory and discuss an equivalence between the two nontrivial vacuum 
solutions in this approximation. In the Appendix B, we present a proof of unitary inequiva- 
lence between the bifurcated nontrivial solutions from the OR method in the infinite volume 
limit. In the Appendix C, we summarize the results of the spurious solutions with some 
discussion. 



II. THE OSCILLATOR REPRESENTATION (OR) METHOD 

The usual recipe to produce a quantum system comes from the correspondence principle. 
Namely, a quantum system can be obtained from its classical analog by assuming that the 
field is not a c-number but an operator satisfying relevant commutation relations. The 
oscillator representation method is a generalization of this scheme. To start with, let us 
consider the 4 theory as a simple example. 

The Hamiltonian density of the classical field is given by 

^ = ^ 2 + i(V0) 2 + im 2 2 + f0 4 . (4) 

To plug the quantum physics in, we postulate that the field (f> and the corresponding conju- 
gate momentum n are operators which satisfy canonical equal time commutation relations 

[<f>(x), n(y)] X0=V0 = i£(x - y). (5) 

The usual approach to quantize the field, i.e. to find the operators 4> and it which satisfy 
the condition (5) is to start from the free field, where g=0. Then the operators are 
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M*)= [ -S= 7 =L= (a m (k)e ik -* + ai(k)e- ik -% 
J V27T v /2cj m (k) 

Ux) = 7 / ^y^KMe** - o+(k) e -**), (6) 

where fc = (cj m ,k), w m = Vk 2 + to 2 and the operators of creation and annihilation of the 
particle of mass to and momentum k satisfy the commutation relations 

[a m (k),a+(k')]=5(k-k'), 
[a m (k), a m (k')] = [a+(k), a+(k')] = 0. (7) 

With this definition the state a+(k)|0) is an eigenstate of both free Hamiltonian and mo- 
mentum operators with corresponding eigenvalues u m and k. 

In the case of interacting fields the states a+(k)|0) are not already eigenstates of the full 
Hamiltonian. Nevertheless, in the case of small interaction strength, the free representation 
of field operators can be a good starting point to develop a perturbation theory. The 
Hamiltonian density operator is then written as 

W-»Wm=^m+ \Wm) 2 + + |<^. (8) 

Now, we use the canonical transformation of the field operators to build the quasiparticle 
representation of the same Hamiltonian. We consider the essential point of the reparame- 
terizations, i.e. the change of mass, which can affect the choice of representation for the free 
field quantization. For the free field operators, the change of mass to — > M of particles is 
described by the Bogoliubov-Valatin transformation: 

if / u> m (k) I vmQs) \ ,,v if j umps) I wmQ*) \ + 

°" (k) = 2 + y ) aJk) ~ 2 [i^m ~ v ) a - { ~ k) - 

+ /, x M /^m(k) uj M (k)\ +f . if u m (k) u M (k)\ . . 

a « (k) = 2 (V^w + V^w J a '» (k) " 2 (V^w " V^w J am( - k) - (9) 

Here, the free field operators with mass M are expressed in terms of the field operators with 
mass to. 

A consequence of the ambiguity of the mass definition due to interactions comes from 
the ambiguity in the choice of the initial representation of the interacting field. This was 
emphasized in Coleman's paper [13]. He also proposed a useful technique on how to redefine 
a normal ordered product of any number of field operators with respect to a new value of 
mass in the case of (1+1) scalar field theory. The redefinition after the mass change is given 
by the formulas [3,13]: 

Nrn Qt£ + i(V0 m ) 2 ) = N M + i(V0 M ) 2 ) + ^(M 2 - TO 2 ), (10) 

N m (e*+") = (—) N M (e^), (11) 
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where (3 is some arbitrary c- number, N m (NM) stands for the normal ordering of an operator 
with respect to the initial mass m (the new mass M), <p m is the free quantum field defined 
by Eq.(6) and 0m is the quantum field of independent quasi-particles for which the repre- 
sentation given by Eq. (6) is valid after changing m — * M. The correspondence between the 
canonical transformation and the renormalization group equation is discussed in Ref. [5]. 
The expression (11) can be rewritten as 

[HI 

/I M 2 \j ( — 1 V 

iuc) = E( s ln <^>) w^m N "^ < 12) 



where [|] is the integer part of 



n 

The most simple illustration of nontrivial solutions in this case is the 1+1 dimensional 
4 theory. Rewriting the Hamiltonian (8) in terms of the quasi-particles using Coleman's 
formula of normal ordering rearrangement, one gets 

N m {H m ) = N m (Hm) 

+ \{m 2 - M 2 )7V M (0 M ) + ^(M 2 - - 2 ) " ^ " |W(<&) + |f£, (13) 

where t = ln(^). 

For a consistent description of the dynamics in terms of the quasi-particles one should 
provide the right form of the Hamiltonian, i.e. the free part should correspond to that 
of the field with mass M and the interacting part should contain only terms with powers 
of (f>M larger than 2 [5]. As it can be seen, the expression (13) itself can not provide this 
requirement for the field with mass other than m. Revising just the mass does not lead to a 
nontrivial solution and thus an additional transformation is obviously required. The latter 
can be done by shifting the quantization point of the field, which makes the quanta of the 
field to be produced around some vacuum expectation value (0) 7^ 0. 

Thus, let us now change the field variable from cf) M (x) to $m(^) for the inclusion of the 
field condensation denoted by b(— (0jw-)):ie., 

Mx) = $ M (x) + b. (14) 

Then the Hamiltonian density is written as [5]: 

N M (H) = N M (HZ ght ) + Hi + e M , (15) 

where TC^ ht is the right form of Hamiltonian density for the field with mass M, 

«f = l«*M + ^(V$ M? + + f <t>M + gb&M, (16) 

Em is the energy density of the vacuum for the quasiparticle of mass M up to the zeroth 
order of the dimensionless effective coupling g/M 2 , 

_ m 2 b 2 gb A m 2 t 3gb 2 t 3gt 2 M 2 - m 2 

£M -~ + T - ¥"17 + 647 + Stt ' (17) 
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and 



Wi = {\{m 2 - M 2 ) + - + {m 2 b + - ^}$ M (18) 

is a remainder. To have the Hamiltonian in the right form the remainder must be equal to 
0. Therefore, we have two equations of two variables, t and b, as follows: 



and 



&(l + Gb 2 



where G is the dimensionless coupling of the initial theory, i.e. G = g/m 2 . 

We note here the correspondence to the GEP discussed in Section I. The energy density 
Em given by Eq.(17) corresponds to < 0o,fi|7Y|fi,0 o > presented in Eq.(2). Minimizing 
< O , 0o > (or s M ) with respect to Vl (or t), one gets the GEP given by Eq.(2) with 

the so-called Cl equation 1 which is identical to Eq.(19) that is one of the OR equations. The 
other OR equation, Eq.(20), can be obtained by making the total (not partial) derivative of 
the GEP (i.e., Eq.(4.11) of Ref. [8] or Eq. (2. 16b) of Ref. [10]) with respect to the condensation 
<f>o(— b) and setting it to zero, i.e. dVa{4>o) / 'd(f) = 0. One should note, however, that the 
GEP is a function of the global field condensate O n °t the local field 4>(x). Thus, the GEP 
is not a part of interaction Hamiltonian describing the field dynamics. In our work, instead 
of using the GEP, we will directly use the interaction Hamiltonian and the corresponding 
classical potential (see e.g. Eq.(22)). 

The OR equations (Eqs.(19) and (20)) have three solutions and one of them is trivial: 
t — 0, b — 0. The two nontrivial mass solutions are shown in Fig.l as solid and dashed 
lines. The dashed line solution on Fig.l becomes similar to the trivial one when the di- 
mensionless coupling becomes large: Mq^oo m, bc^oo 0. However, the solid line 
solution from Fig.l is what we are interested in because it fulfills the duality. For this case, 

Mq^^ — > ^m 2 G\nG, &g->oo \J In G. As one can see in this solution, when the di- 
mensionless coupling G of the initial theory becomes large, the dimensionless couplings of 
the quasiparticle theory, — gb/M 2 and — g/M 2 , get small [5]: XgLoc ~^ V^g 1 
XgLoo ~ *■ 3Gi g ■ Thus, the quasiparticle theory becomes perturbative when the initial cou- 
phng G gets very strong. The physical mass in this regime is close to M « and 

the corrections are in the order of 0((x^) 2 )- I n this way, we can see the duality between the 
initial nonperturbative theory of unbroken — > —0 symmetry with a large dimensionless 
coupling G and the quasiparticle theory of broken $ — > — $ symmetry with small couplings 



x See Eq.(4.12) of Ref. [8] which was rederived as Eq.(2.16a) of Ref. [10] using Bogoliubov 
transformations. 
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and ■ This duality is exact because the canonical transformation leaves the Hamil- 
tonian intact but changes only its representation in terms of the quasiparticle field variable 
$ M (x). 

As we discussed earlier, the multiple solutions are not unexpected because the OR equa- 
tions, Eqs.(19) and (20), are algebraically nonlinear. However, it can be shown (see Appendix 
B) that the two vacua associated with different values of t or b are unitarily inequivalent in 
the infinite volume limit. Thus, the two nontrivial solutions are not connected to each other. 
The bifurcation of quasi-particle masses appear above the coupling G ~ 9.04. However, the 
correct critical coupling of the phase transition should be determined from the comparison of 
the energy density between the trivial and condensed vacua. In Fig. 2, we present the energy 
density of the nontrivial solutions and find that the condensed vacuum energy density of 
the duality-related solution gets lower than the trivial vacuum energy density only above 
G ~ 10.21. This value G ~ 10.21 coincides with the critical coupling constant obtained by 
the GEP method since our energy density given by Eq.(17) is equal to the GEP at the min- 
imum. This value is remarkably close to the critical coupling constant (10.24(3)) obtained 
by the lattice calculation [14]. Nevertheless, we note that our OR result doesn't take into 
account the the higher order corrections in the quasiparticle effective coupling g/M 2 and 
obtain the first order phase transition instead of the second order phase transition which 
may be more accurate description of the phase transition in (0 4 ) 1+1 theory [11,15]. In this 
work, we thus focus only on the region where the original coupling is much larger than the 
critical coupling so that the perturbative calculation in terms of quasiparticle degrees of 
freedom is allowed. 

Although there is a coincidence of the critical coupling, it is not clear if the values of 
quasi-particle masses defined in the GEP, HA and OR methods also coincide. In the GEP 
approach, it seems only natural to define the quasiparticle mass as the second derivative 
of the GEP 2 , i.e. g? 2 Vg(0o) / dcf^. As discussed earlier, in the GEP approach a variational 
method is taken to minimize the GEP satisfying just Eq.(19). Confirming the well-known 
equivalence between the GEP method and the HA for the ground-state energy [1,6], the 
GEP result using this definition of quasiparticle mass, i.e. the second derivative of the GEP, 
indeed agrees precisely with the Hartree results (squares and circles in Fig.l) 3 . However, 
one should note that the second derivative of the GEP at the minimum O = b does not 

yield the mass M of the quasi-particle precisely but somewhat modified value Myj ^jp^§|; 

i.e. with the notations used in Eq.(4.11) of Ref. [8], the second derivative of the GEP, V G , 
is given by 



2 See for example Eq.(3.2) of Ref. [8]. 

3 In Chang's paper [11] he found two different types of nontrivial solutions using Hartree method 
(47rc 2 > 3 and 47rc 2 < 3 in his notations). Nevertheless, the two type of solutions lead in fact to 
the same physical results as shown in Fig.l. For more details, see the Appendix A. 




(21) 
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where x = A# = ^2 and O — b in our notations. Removing t(— lax) in Eqs.(19) and 
(20), one can immediately get b = identical to O = ±^/^|-. This is the reason why 

the GEP result given by Eq.(21) shown by dotted line in Fig.l doesn't coincide either of 
the bifurcated OR solutions (solid and dashed lines). Of course, if both the fl equation and 
the GEP minimum condition dVa((f)o) / d<fio = are solved simultaneously to determine f2, 
then it is identical to the OR method and the value of fl coincides with the quasiparticle 
mass M defined in the OR method. However, we note once again that the GEP and HA 
methods are quite different from the OR procedure because the GEP and HA methods rely 
on the variational procedure while the OR method is explicitly solving both Eqs.(19) and 
(20) simultaneously. 

As we have already discussed, the solid line in Fig.l exhibits a duality between the theory 
with mass m based on the trivial vacuum 6 = and the theory with mass M based on the 
nontrivial vacuum b 7^ 0. Namely, a strong coupling theory with mass m is identical to 
a weak coupling theory with mass M. A similar observation of duality can be found in 
the solutions obtained by GEP and HA (see also Appendix A), although the quasiparticle 
masses defined in the GEP and HA approaches do not coincide with the quasiparticle mass 
defined in the OR method. We also note that the GEP result (dotted line) approaches to 
the OR result given by the solid line in the limit G — > 00 since the terms of x overwhelm 
the terms of As in Eq.(21) for that limit. 

In Ref. [10], it was argued that the method of Bogoliubov-Valatin transformations gives 
the same results as the GEP. In fact all of the nonperturbative methods that we discuss in 
this work,z.e. OR, GEP and HA, can be described formally in terms of Bogoliubov-Valatin 
transformations because all of them use only the change of particle mass and the field shift 
to get the new representation of the field variables. However, these transformations contain 
parameters (mass ratio and the value of the shift) which can not be determined from the 
transformations themselves. To fix the parameters, some additional requirement is necessary. 
The authors of Ref. [10] used essentially the same variational requirement as the GEP and 
therefore it is very natural that they found the same results as the GEP. On the other 
hand, in the OR method, we are not using the GEP at all but relying on the form of the 
Hamiltonian described in Section I. 

Applying the correspondence principle to the Hamiltonian (16), we now consider the 
classical potential given by 

yclass^class = ^class _ ft) = _^class _ ft) 2 + 9_ {(j) class _ fe) 4 + potass _ &)3> (22) 

which is written here in terms of the classical field corresponding to the initial, unshifted 
field 4>m- As it can be easily seen, the potential has extremum at (f) c ^" ss = b. This is provided 
by the demand not to have terms in Hamiltonian which are linear in $m- The double 
derivative in this point is larger than and equal to M 2 , i.e. the classical potential has 
minimum at <p C M SS = b and the mass corresponding to oscillations around the minimum is 
precisely M. That is why the OR method provides self-consistent procedure of quantization 
of interacting fields and interpretation of the solutions as particle excitations is valid. (In 
contrast, ^j^-\s =b 7^ f° r the GEP as shown in Eq.(21).) In addition, one can easily 
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prove the symmetry of V dass under the exchange of <p C M SS and —(p c ^ LSS ,i.e.V class ((l) c ^ ,ss ) = 
ydass^dass^ ut iii z i ng the OR equations (Eqs.(19) and (20)). 

In Fig. 3, the effective potential V class given by Eq.(22) is presented. As one can see in 
Fig. 3, the bifurcated solutions for G > 9.04 indeed generate two effective potentials (I, II) at 
a given coupling G, e.g. G = 10 or 12. Both potentials (I and II) in Fig.3 have clear minima 
and the second derivative at each minimum generates the corresponding quasi-particle mass 
M precisely. The growth of bifurcation is also evident from the comparison of the potentials 
I and II for G — 10 and 12. We note however that only the potential I provides the physically 
meaningful duality-related solution and the potential II doesn't satisfy the duality that we 
discuss in this work. 

An analysis of realistic models using the oscillator representation method is more com- 
plicate (see, for example, discussion for 3+1 dimensional 4 theory in [5]) due to substantial 
sophistication of the renormalization procedure. Even the simple models however were not 
fully studied in the framework of OR method. As presented in this Section, the comparison 
of the OR predictions with the GEP and Hartree results has been done for 1+1 dimen- 
sional 4 theory [7,11]. In higher dimensions, the variational methods such as the GEP 
and HA may have a rather fundamental difficulty The reason of this problem originates 
from the error estimates of the approximation in the GEP and Hartree methods. For in- 
stance, the HA replaces the interaction term A0 n for n > 2 by the kinetic term Am 2 2 + E . 
The task of getting the better approximation is to minimize the positively defined mea- 
sure ||A0 ra — Am 2 2 — Eq\\. Nevertheless, the vacuum expectation value of the operator 
\(f) n — Am 2 (p 2 — Eq is not positively defined. Thus, it is necessary to calculate the vacuum 
expectation value for its square that is always positive. In that case, even if the interaction 
term itself is renormalizable, its square may be not (in particular, for D > 2). That is why 
in higher dimensions the GEP and Hartree methods come across a renormalization problem 
in computing their error measures [16]. More discussion about the difficulty in the GEP 
method for the higher dimension can be found in Ref. [17]. However, this is not the case for 
the OR method, since it does not use any variational procedure. 

Therefore, among the three (OR,GEP,Hartree), the OR method seems to be most at- 
tractive. The reasons may be summarized as follows: 

i) The problems of renormalization in higher dimensions are technical, but not principal. 

ii) The Hamiltonian is identically rewritten via canonical transformation and consequently 
the transformed Hamiltonian keeps the identical information as the original Hamiltonian. 

iii) The second derivative of effective potential gives the quasi-particle mass precisely in the 
OR, while the others give somewhat modified values. 

III. APPLICATION OF THE OR METHOD TO THE {4> 6 )i+i THEORY 

Using the OR method, we now compute the parameters of different phases for the scalar 
field theory with both 4 and c/> 6 interactions. The starting classical Hamiltonian density is 

n m = \n 2 m + i(V0 m ) 2 + ^ + 9 -4>t + Hi- (23) 

The way to proceed is the same as in the case of 4 theory. Using the formulas (10) and 
(11) and shifting the quantization point, we get: 
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N m (H m ) = N M (HZ 9ht + U x ) + e M , (24) 
where the "right" Hamiltonian density 7i r ^ ht is given by 

nT = I^m + ^(V$m) 2 + ^ 2 M + 9 ~y M + + gb** M (25) 



2 2 V ' 2 4 



15/1*,. 15/ifa 



+ 6/i&$^ + 15/ifo 2 ^ + 20/i6 3 $^ - — $ 4 M 

47T 7T 

the remainder TL\ is given by 

( l , 2 3^6 2 irl4 , 456 2 /it 45/it 2 \ 2 

= -(m 2 - M ) + — — + lhb 4 h + $L 

V2 V ; 2 8tt 2tt 16tt 2 / m 

/ 2i ,, 3gbt - 15/i6 3 t 456/it 2 \ . . 

+ (m% + gt?-JL- + 6hb 5 - — — + -^-) $ M , (26) 

and the energy density Em of the vacuum is given by 

m 2 b 2 gb 4 m 2 t 3gb 2 t 3gt 2 



£ M = — ; V — h 



2 4 8tt 8tt 64tt 2 

1 , - 2N , j6 15/i6 4 t 45/i6 2 t 2 15/it 3 
+ _ (M ._ m » )+ ^ __ + ___. (27) 

As in the last section, here we use the notation $ M = 4> M — b. 

The requirement to have the Hamiltonian in the right form means that the remainder 
Hi is equal to zero. This gives the following two equations for t and B (B = b 2 ): 

, 3Gt r „ AhHBt 45Ht 2 . . 

1 - e* + 3G5 + 30iLB 2 + — = 0, 28 

47T n 8n 2 



3Gt 9 IbHBt A5Ht 2 , x 

1 + GB - — + 6iLB 2 + — — = 0, 29 

47T n 8n 2 

where G = g/m 2 and H = h/m 2 are the dimensionless couplings. The conditions (28) and 
(29) again provide the minimum of the classical potential at 4>f^" ss = b and the right value 
of the double derivative around it. 

Solving Eqs. (28) and (29) numerically, we found five nontrivial mass solutions shown 
in Figs. 4, 5 and 6. The corresponding vacuum energy densities are presented in Figs. 7, 
8 and 9. In Fig. 7, two identical plots are shown in two different angles of view to make 
clear presentation of the result. Among the five, two solutions correspond to the symmetric 
phase, i.e. b = 0. They are represented by the upper and lower branches (or sheets) of 
the manifold shown in Fig. 4. The vacuum energy densities for these solutions are shown in 
Fig. 7. The upper and lower branches in Fig.4 correspond to the lower and upper branches in 
Fig. 7, respectively. The upper branch in Fig.4 is the physically meaningful duality- related 
solution while the lower branch in Fig.4 does not provide the duality that we are interested 
in this work. The duality-related solution has the lower energy density as shown in Fig. 7. 
The dimensionless couplings which correspond to the <J> 6 and $ 4 interactions of the duality- 
related quasiparticle theory are given by 
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x w = A x (4) = _9_ _ m 



Note here that x^ includes both couplings of g and h. This is due to the self-organization of 
fields via the quantum fluctuations in the vacuum, i. e. the original 6 interaction induces the 
effective $ 4 interaction with the negative effective coupling in this case. Their dependencies 
on H are presented in Fig. 10, where the original 4 interaction coupling g is put to zero for 
simplicity. It can be seen from Fig. 10, both x^ and x^ become small as H grows. Thus, 
the qusiparticle theory in the domain of large H can be solved by a standard perturbation 
method. 

Among the rest three solutions, the two solutions corresponding to the broken-symmetry 
phase with the nonzero real condensation {b = \^B) are shown as the two branches in 
Fig. 5. The vacuum energy densities for these solutions are shown in Fig.8. Again the upper 
(lower) branch of Fig. 5 corresponds to the lower (upper) branch of Fig.8. Similar to the case 
of symmetric phase (Fig.4), the lower branch in Fig. 5 does not provide the duality that we 
are interested in this work. However, the upper branch in Fig. 5 corresponds to the duality- 
related solution and gives the lower value of the vacuum energy density (the lower branch 
in Fig.8). For this broken-symmetry phase, the dimensionless couplings of the quasiparticle 
interactions are given by 

( 6 ) h (5) 6hb ( 4 ) g 15hb 2 15ht 

(3) _ gb 20hb 3 15hbt 
Xbs Ml + Ml nMl ■ [6 } 

The if-dependence of these couplings is shown in Fig. 11 (again, we put g — to simplify the 
presentation). At sufficiently large H the couplings get small and the quasiparticle theory 
becomes perturbative. When the dimensionless couplings of the symmetric solution go to 
zero, X^s\x^ — ► 0, the dimensionless couplings of the broken-symmetry theory also go to 
zero as shown in Fig. 12. 

The last solution shown in Fig. 6 has a pure imaginary value for the condensate b {i.e. 
negative B). In Fig. 9, we show the corresponding vacuum energy density. Because the 
Hamiltonian in this case becomes non-hermitian, we call these solutions spurious. The spu- 
rious solutions appear only after the inclusion of 6 interaction. Thus, they are consequences 
of the higher nonlinearity in the OR equations. Since they don't seem to bear any interesting 
duality property that we discuss in this work, we don't present their results in the main text 
but just summarize those in the Appendix C with some discussion. 

The most interesting feature of the 6 theory is the appearence of the two apparently 
different quasiparticle perturbation theories (see the upper branches of Figs. 4 and 5) which 
represent the initial nonperturbative theory. The first one corresponds to the symmet- 
ric phase with the dimensionless quasiparticle couplings given by Eq.(30) and the second 
one corresponds to the broken-symmetry phase with the dimensionless couplings given by 
Eq.(31). The subscripts "s" and "bs" of the coupling x an d the mass M represent the corre- 
spondence to the symmetric and broken-symmetry phases, respectively. The two quasipar- 
ticle theories {i.e. Xs's and Xt>s' s ) differ from each other significantly in the sence that one of 
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them (x s 's) preserves the — > —0 symmetry of the initial Hamiltonian (23) upon changing 
$ — > — $, while the other (x& s 's) does not respect the $ — > — $ symmetry. Physically, this 
means that in the theory of symmetric phase with the Xs couplings the processes involving 
odd number of particles are not allowed while they are allowed in the broken-symmetry the- 
ory with Xbs couplings. This is because in the theory of symmetric phase, the Hamiltonian 
depends on the even powers of <£> (to provide $ — > — $ symmetry) and thus the S-matrix 
contains only the even number of creation and annihilation operators. 

However, one may wonder if these two apparently different quasiparticle theories are 
in fact equivalent to each other describing the same physics with just different degrees of 
freedom. In this respect, we notice that the two different quasiparticle theories give the 
same effective potential in the limit of very large G and H. In Figs. 13, 14 and 15 we show 
the classical potentials which correspond to the two quasiparticle theories at different values 
of the coupling x^ ■ As it can be seen, when the quasiparticle coupling x^ becomes small 
(i.e. H becomes very large), the relative difference between the two quasiparticle potentials 
diminishes substantially. We verified that the two classical potentials are indeed identical 
to each other in the x^ ~^ limit. 

In the GEP approach [9], it was also recognized that there are three solutions of sym- 
metric phase (one of them is the trivial solution and two others are the nontrivial solutions 
corresponding to the ones that we showed in Fig. 4). As we consider only one of the two 
nontrivial solutions (the upper branch in Fig.4) as the physically relevant (duality-related) 
solution, the GEP approach also dealt with only one (physically relevant) solution among 
the two nontrivial solutions. In Ref. [9] it was shown that the form of Q equation of the 
GEP approach is invariant under the exchange of m and M s with an appropriate redefinition 
of effective couplings G' — (G — ^-t )e~ to and H' = He~ to , as well as t' = t — t , where 

t — In ^f- satisfies the constraint equation given by Eq.(29): 1 — e~ to = — 4 g^° (Note 
B = for the symmetric phase). Indeed, we can show the form invariance of both OR Eqs. 
(28) and (29) under the same transformations. 

Observing such form invariance, it was argued in the GEP approach [9] that the nontrivial 
symmetric phase solution is just a duplication of the trivial solution and thus should be 
avoided. However, as we discussed above, the duality between the trivial and nontrivial 
symmetric phase solutions provides the physical meaning to the form invariance of the OR 
equations. Because of the duality, we see the utility of both trivial and nontrivial solutions in 
the symmetric phase. While the trivial solution can be utilized for the perturbation theory 
at small G and H, the nontrivial symmetric phase solution can be used for the construction 
of quasiparticle perturbation theory in the regime of very large G and H. 

Furthemore, since there exists also a duality between the trivial solution and the nontriv- 
ial broken-symmetry solution we can rather easily understand the equivalence of the classical 
effective potentials between the symmetric phase and the broken-symmetry phase for the 
very large couplings of the original theory (see Fig. 15). In the limit of x^ ~ * (which 
means also Xbs's-^ 0, as shown in Fig. 12), the higher order corrections get negligible and 
thus the equivalence of the classical potentials between the symmetric and broken-symmetry 
phases is revealed manifestly as shown in Fig. 15. 
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IV. SUMMARY AND CONCLUSIONS 



The formation of quasiparticle in RQFT may be a good example of self-organizing nature 
in the nonperturbative strong interactions. This nature is realized by altering the vacuum 
structure and lowering the vacuum energy density as the strength of the interactions gets 
more enhanced. Our main goal in this paper was to provide the explicit examples that 
reveal such nature of self-organizing RQFT and show the utility of duality generated by the 
canonical transformation in solving the nonperturbative strong interaction problem. For 
the explicit but simple examples, we analized of 1+1 dimensional 4 and c/> 6 theories. Our 
purpose was fulfilled by applying the OR method which is a generalization of the canonical 
quantization scheme. 

The OR method uses canonical commutation relations for the free field operators and 
employs the canonical transformation of quantum fields represented by the change of the 
field mass and the condensation of the field. The idea to demand the right form of the 
Hamiltonian, i.e. excluding the terms linear in $ and taking the coefficient of $ 2 as M 2 /2 
guarantees that the classical potential has the minimum at the quantization point 4>f^" ss = b 
and that the double derivative of the classical potential at this point is given by M 2 . This 
assures that the OR method gives a self-consistent procedure of quantization and a valid 
description of the field excitations around the minimum point in terms of quasiparticles. 

In contrast to the GEP and HA methods, the OR method is not based on the varia- 
tional procedure so that it does not have a bottleneck of renormalization problem in higher 
dimensions as discussed in Section II. However, the OR method has not been utilized as 
extensively as the GEP and HA methods. Even in the 1+1 dimensional scalar field theory, 
not much of the OR results were contrasted from the GEP and HA results although the OR 
results for the (0 4 )i+i theory were presented in Ref. [5]. We attempted to fill such gap in 
this work. 

For the fixed classical Hamiltonian and renormalization procedure (fixed by the counter- 
terms), it is natural to expect that the OR method yields several different solutions for 
the quasiparticle mass and condensation because the OR equations are in general nonlin- 
ear. However, the unitary inequivalence arising from these different OR vacuum solutions 
prohibits the transitions among the quasiparticles with different mass and/or condensation. 
Although some of the solutions correspond to the heavier quasiparticles, they do not decay 
into the light ones due to the unitary inequivalence; i.e. each sort of the quasiparticles is 
stable within the OR. The unitary inequivalence between the OR solutions is shown in Ap- 
pendix B. A recent application of the unitary inequivalence to the flavor mixing phenomena 
can be found in Ref. [18]. In this work, we selected only the duality-related OR solu- 
tions because they can provide the immediate physical interpretation. The reason why the 
physical interpretation is immediate for the duality-related solutions is because the strong 
coupling particle theory becomes solvable once it is transformed into an exactly equivalent 
quasiparticle theory with weak couplings. 

Using the (0 4 )i + i theory, we compared the OR results with the results obtained by GEP 
and HA. In Ref. [5], it was discussed that the value of the critical coupling G ~ 10.21 
obtained by the OR method coincides with the one obtained by the GEP method. We also 
note that this value of the critical coupling agrees remarkably well with the lattice result 
[14] . Due to the equivalence between the HA and GEP methods for the ground state energy 
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[1,6] the agreement of the critical couplings between the OR and HA methods follows as 
well. Although the critical values of coupling constant G turned out to be same for all three 
methods (OR, GEP, HA), the OR method gives the quantitatively different dependence of 
the nontrivial field mass on the effective coupling constant, compare to the GEP and HA 
results (see Fig.l and the discussions in Section II). 

We have extended the application to the 1+1 dimensional 6 theory and compared 
the results of OR method with what was already obtained by the GEP. We found that 
the different solutions for the mass of quasi-particles between the OR and the GEP are 
sustained. In the OR analysis of (0 6 )i + i theory we found two physically meaningful duality- 
related solutions (upper sheets of Figs. 4 and 5), one in the symmetric phase and the other 
in the broken-symmetry phase. Although these two solutions have apparently different 
symmetry property, we find that their effective potentials do agree in the limit of very large 
couplings G and H of the initial theory, i.e. very small effective couplings Xs' s an d Xbs' s 
of the quasiparticle theories. This indicates an existence of a unique effective potential to 
the quasiparticles regardless of their symmetry properties. The actual perturbation theory 
using the quasiparticles degrees of freedom with the systematic higher order corrections may 
deserve further investigation. 
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APPENDIX A: Hartree Approximations in the (0 4 )i+i Theory 

The essential idea of the Hartree approximation [11] is to linearize the non-linear field equa- 
tions using the mean fields, e.g. 3 — > 3(0 2 )0 — 2(0) 3 7 etc., splitting the fields into classical 
mean fields (0 C = (0)) plus the quantum fields(0 g ), i.e. = C + <p q . In the classical limit, 
we know that the ground state is given by 2 = c 2 if c 2 is defined by c 2 = — ^- and pos- 



itive (i.e. c > 0). The quantum fluctuations {(jr q (x)) are found to be infinite. However, 
one can renormalize the theory using a mass counter term \B(p 2 1 where the constant B 
is determined by requiring that C = c is a static solution (corresponding to the so-called 
"abnormal vacuum state" or a sort of vacuum condensate) of a Hartree equation and this 
leads to B = — 3g(0 2 )^ c=c . After renormalization, the effective potential V(4> c ) is obtained in 
terms of the subtracted quantum fluctuation, i.e. A(0 2 ) = (0 2 )0 C — (0 2 )0 C = C , which satisfies 
the iterative (mass-gap type) equation [11]; 

A{ ^-i ln ^-f +3 A^y (A1) 

Using the solutions of Eq.(Al), one can find that <p c = c is a local minimum (maximum) of 
V(<f> c ) for 47rc 2 > 3(0 < 4ttc 2 < 3) while C = is always a local minimum. For < 47rc 2 < 3, 
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a nontrivial local minimum exists at <f) c = ^jc 2 — 3A(0 2 ) — c ' . Expanding V((f) c ) in each of 
the local minima, one can introduce the following mass parameters; 



2 d?V(<j> c 
mi - 



°" # 2 



o = ^(3A(0;)U c= o-c 2 ), 



2 _ d 2 V(<Pc 



2gc 



8ttc 2 - 6 
8ttc 2 + 3' 



m 2 



d2V ^c) A , 2 oA/^^ 47r (c 2 -3A(0 2 ))-3 



# c 2 



4^-3A^)) Mc2 _ 3A -- +3 . (A2) 



Then, the intrinsic strength measured in terms of these mass parameters are given by the 
dimensionless couplings; i.e. 

9 9 9 , AQ x 

9o = 9c = — 9c' = (A3) 

wig m l c m z c . 

For the entire c 2 > region, it turns out that the nontrivial vacuum solutions exist only 
f° r <7o > fi'crit ~ 9.045 (note that it is not the same as the critical coupling of the phase 
transition G = 10.21 as discussed in Section II), where the critical coupling g = g crit occurs 
at 47rc 2 = 3. Thus, it appears that one can find the two nontrivial vacuum solutions; one 
for 47r c 2 > 3 and the other for < 47rc 2 < 3. However, we find that the two solutions are 
indeed identical in the sense that g c and g c > coincide when we plot them (see Fig. 16) as a 
function of g , i.e. g c (g ) = g c >(g ) for all g > g crit . 

We also notice a duality between the theory with g based on the trivial vacuum (f> c = 
and the theory with g c (or g c ,) based on the nontrivial vacuum C = c for 47rc 2 > 3 (or C = d 
for < 47TC 2 < 3). Namely, a strong coupling g theory is identical to a weak coupling g c 
(or g c t) theory. 

APPENDIX B: Unitary Inequivalence between the Nontrivial OR Solutions 

In this Appendix, we show the unitary inequivalence between a vacuum with b\ and t\ 
and another vacuum b 2 and t 2 , where t-y ^ t 2 or b\ ^b 2 . 

The canonical transformation taking a trivial vacuum {i.e. b = and t — 1) into a 
nontrivial one {i.e. b ^ and t ^ 1) has the following form {e.g. for the bosons) 

1. afc — > ak — 2nmb5{k); 

2 ak — > ak cosh(A) — a_ k sinh(A), (Bl) 
a k — > a k cosh(A) — a_fc sinh(A). 

Their corresponding operators have the form 

Ui = exp < —2iimb{a — a ) > ; 

r - t t l ( B2 ) 

U 2 = exp <±f dk\{k){a_ k ak - a\ a -k) f 
with \{k) = \ In ^gy and w(ife) = VA: 2 + t 2 m 2 . 
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To connect the two vacua, one with b\ and t\ and another with b 2 and t 2 , we can 
use the transformation given by Eq.(B2) with b = b 2 — 61 and \(k) = | In w here 

uJi{k) = \Jh 2 + t 2 m 2 . To show that, let us define V; L : (0, 1) — > {h,ti) by the transformation 
given by Eq.(B2). Then 

(Mi> -> (M2) = (&i,ti> -> (0,1) -> (M2) 

and therefore F = V 2 Vf 1 : (61, ii) — > (b 2 ,t 2 ). Straightforward algebraic manipulation with 
Eq.(B2) then results in the unitary inequivalence between the vacuum with (6 1? t x ) and the 
vacuum with (b 2 ,t 2 ). 

In Ref. [5], it was mentioned that 

c/ 2 = e -(^/^ ln ( cosh [ A ( fe )]) . ^ . 

where [a , aj] = <^(0)- Here, we can equate <5^ 3 ^(0) to t^w- Note also that transformation 
C/i acts only on the k — field operators, therefore 



t/ 2 o[/i= n^) Mo^o)), 

( II ■U 2 (k):) : C/ 2 (0) o ^(0) : . 



: C/ 2 o 
Hence 

U 2 oU 1= e ( ' fj : U 2 (k) : (C/ 2 (0) o U^O)) . 

k^O 

Then for the inner product of the vacuum |0) 1 associated with (bi,ti) and |0) 2 associated 
with (b 2 , t 2 ) we have 

1 (0|0) 2 = n i,fc (0|0) 2 , fe < II U (°l°>2,fc = e ( ' ^ - 0, 

k^O 

as 

V -> 00. 

The two vacua therefore are unitary inequivalent if their t parameters are different. If t\ = t 2 
but 61 ^ b 2 then \(k) = and V = C/i(6 2 - 61). Then 

1 (o|o) 2 = J] u (o|o) 2 , fc = i,*=o (o|o) 2;fc=0 = e-^?m = 0, 

in the infinite volume limit. 

Therefore 1 (0|0) 2 = in the infinite volume limit whenever t\ 7^ t 2 or b\^b 2 . 

APPENDIX C: Spurious Solutions in the (0 6 )i+i Theory 



17 



Among the solutions of OR Eqs.(28) and (29) for 1+1 dimensional 6 theory there is one with 
negative B that gives pure imaginary condensate b = i\J\B\. This solution is shown in Fig. 6 
and the corresponding vacuum energy density is shown in Fig. 9. In this case the Hamiltonian 
becomes non-hermitian because it contains anti-hermitian terms in odd powers of <3>m- When 
one calculates average values of such operators in any state with definite number of particles, 
these terms drop out, because odd number of creation and annihilation operators can not 
leave the number of particles unchanged. However, the physical states are complicate and 
contain indefinite number of particles so that the antihermitian terms would contribute to 
the physical observables and give them imaginary expectation values. Moreover, our starting 
Hamiltonian density given by Eq.(23) involves only the real (not complex as in the 0(2) 
symmetric theory) scalar field <p m and even after shifting the quantization point the field 
$m is required to remain as real not complex. For these reasons, we call the imaginary b 
solution spurious. In principle, we can throw away the spurious solution and do not consider 
the imaginary-?) phase in the simple one-component model. 

Nevertheless, in the history of many-body analysis sometimes the occurance of such spu- 
rious solution (forbidden in principle) led to interesting physics. For example, the occurance 
of imaginary solution has been previously observed in the analysis of many-body system 
when the method gets more accurate but the number of degrees of freedom is not sufficient. 
Such example may be illustrated in the random phase approximation (RPA) calculation of 
a fermion system [19] where the imaginary energy poles are found if the fermion interac- 
tion is attractive and weak below a certain critical temperature. These imaginary poles 
were interpreted [19] as the mathematical manifestation of the instability which leads to the 
superconducting state. Such imaginary poles were however not found in the less accurate 
method such as the Tamm-Dancoff approximation (TDA). Thus, if there is any possibility 
that the imaginary solution could actually mean something interesting in the many-body 
analysis, we should not overlook such a possibility. 

Due to such a tantalizing possibility that the imaginary solution could occur indicating 
an existence of a new phase, we note a few more particular features of the spurious solutions 
instead of simply throwing those solutions away. In particular, it may be interesting to note 
that the corresponding quasi-particle mass values become quite small (M « m) when the 
imaginary b solutions have the least vacuum energy density as shown in Figs. 6 and 9. We 
do not know yet, however, whether these spurious solutions are anything to do with the 
condensation of (zero-mass) Goldstone-bosons if the degrees of freedom are added in such a 
way that the fields become complex and a continuous symmetry (see however the discussion 
for 0(2) symmetry near the end of this Appendix) is restored. With the limited degrees of 
freedom at hands, the theory we consider in the present work does not lead to an appearance 
of some effective charge in the theory. For the spurious solutions, the point of quantization 
is fixed as pure imaginary and can not be rotated by any arbitrary phase. Also, in the 
effective Hamiltonian given by Eq.(25), b can be thought as a part of effective coupling. For 
example, consider the term gb§ 3 M in Eq.(25) and build up an effective potential between 
two § M particles intermediating another <& M as an exchange-boson. Then, the coefficient gb 
of § 3 M is an effective coupling of the triple vertices and the square of the coefficient would 
determine an overall sign of effective potential. The detailed form of the potential is of course 
determined only after taking into account the propagator of exchange-boson. However, the 
overall sign of effective potential determines whether the force between the two <& M particles 
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is attractive or repulsive. Due to the sign change of the potential, if some bound states may 
not be formed (although they could be formed for the real b), then the corresponding energy 
eigenvalues would appear as complex values indicating the instability of those states. In 
some sense, the imaginary b resembles a sort of optical potential which has been used rather 
frequently in the phenomenology of many-body nuclear physics. A recent work showing an 
example of how the imaginary coupling can be used in practice was presented in Ref. [20] . 
To show an overall picture of the effective potential in (0 6 )i + i theory, we present the 

dimensionless effective potential V(<j)) = V class (<j)f!j iss ) +Em /m 2 for the broken-symmetry 

solution and the spurious solution with H/G = — 0.209(o; = —f3 in notations of Ref. [9]) in 
Fig. 17 for several different G values. For simplicity, we didn't plot V((f>) for the symmetric 
phase (b = 0) solution. In Fig. 17, denotes 0^ ass and the symmetry V class (4>f^ ss ) = 
V class {-(t>ff ss ) is manifest for the OR solutions with B = b 2 > 0, i.e. Im(<f>) = 0. Indeed, 
no matter what the sign of B is, we can in general prove the symmetry V(4>) = V (—(/>) 
using the OR equations given by Eqs.(19) and (20). Thus, for the spurious solutions drawn 
off from Im((f>) = plane, one can realize the symmetry V((f>) = V (—(f)) by considering 
both planes of Im{4>) — ±|6| although only the Im{4>) = +\b\ solutions are shown in Fig. 17. 
The solid, long-dashed, short-dashed, dotted and dot-dashed lines stand for G—-0. 838,-1. 68,- 
2.51,-3.35 and -4.19 (a=-0. 1,-0.2,-0.3,-0.4 and -0.5 in notations of Ref. [9]), respectively. For 
the coupling constants G = -2.51 (short-dashed), -3.35 (dotted) and -4.19 (dot-dashed), 
the effective potentials corresponding to the three nontrivial solutions including a spurious 
solution with the imaginary b are shown in Fig. 17: e.g. three dot-dashed lines (the spurious 
one is on the plane Im((fi) ^ and the other two real broken-symmetry solutions are on the 
plane Jm(0) = 0) are shown for the single value of G — —4.19. For the comparison with 
the trivial solutions, we added the solid line for G = —0.838 and the long-dashed line for 
G = —1.68 which have the minimum only at Re(<f>) = Im(<f)) = 0. 

In Figs. 18-20, the effective potentials for the spurious solutions are shown. Since the 
spurious solutions have Im((f)) ^ 0, the corresponding effective potentials V(<f>) are complex. 
The real parts of effective potentials for the spurious solutions of G =-0.838(solid line), 
-1.68 (long-dashed line), -2.51(short-dashed line), -3.35(dotted line) and -4. 19 (dot-dashed 
line) are shown in Figs. 18 and 20 while the imaginary parts for the same coupling constants 
are shown in Fig. 19. Also, the horizontal axes of Figs. 18-20 are given by the real field 
Qdass _ (j)dass _ ft ra ther than the complex field 4> c ^ ss due to the imaginary b. However, 
the effective potentials for G = -2.51, -3.35 and -4.19 in Figs. 18-20 correspond to the ones 
drawn off from the plane Im(4>) = in Fig. 17. As shown in Fig.9, the energy density of the 
spurious solution gets the lowest among all the OR solutions for the small coupling such as 
G = -0.838 and -1.68. Especially, the quasi-particle mass gets close to zero as G becomes 
close to zero. This behavior may be noticed from Fig. 18 where the second derivative of 
Re{V class ) near <p c ^ ss = b is the smallest for G = -0.838. The corresponding behaviors of 
Im(V class ) are shown in Fig. 19. It is manifest from Fig. 19 that Im(V class ) contains only the 
odd powers of the real field $^f ss (= <p c ^ ss — b), which can also be explicitly shown from our 
effective Hamiltonian given by Eq.(25). The large field behaviors of Re(V class ) are shown 
in Fig. 20. We can see that the minima of Re(V dass ) appear in the large real field ^ c ^ ss for 
the spurious solutions, which one can also notice from the effective Hamiltonian (Eq.(25)) 
by keeping only the even powers of $m- 
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In Fig. 21, we show the domain of nontrivial solutions that have the vacuum energy- 
density less than zero in the coupling space of G and H . Without the spurious solutions, 
the domain of nontrivial OR solutions coincides with the GEP result (Fig.l in Ref. [9]) 4 
denoted by the solid line of boundary. Taking into account the spurious solutions, however, 
the domain of trivial solutions reduces to the shaded region of the phase diagram shown 
in Fig. 21. The appearance of the spurious phase (B = b 2 < 0) in the domain where 
couplings are small may be interesting because it is quite opposite to the appearance of 
the normal nontrivial phase (B > 0), for which one would expect that the interaction term 
that leads to the formation of vacuum condensate is dominant at the larger coupling and 
therefore the nontrivial phases appear in the domain of large couplings. The presence of the 
nontrivial solution for any values of the couplings is also interesting in connection with the 
Haag's theorem [21], according to which, if the vacuum state of a Lorentz invariant theory 
is unique, then all of its observables are identical with those of the free theory. Even very 
small admixture of 6 interaction to the 4 theory allows to have a set of vacuum states for 
any value of the <f> 4 coupling and thus to overcome the conditions of the theorem. 

In order to explore whether the spurious solution can be considered as an evidence of 
some real solution by adding more degrees of freedom, we applied the OR method to the 
two-component (0 6 )i + i theory. In this model, the imaginary condensation of one component 
could be represeneted by a real condensation of the other component. However, we found 
that the imaginary-6 solution of the one-component model corresponds to imaginary-^ solu- 
tion of two-component 0(2) model, but not to the real solution. The matching was done by 
solving an arbitrary two-component model 

/ 1 I th 2 

n = iv m (-K 2 + 4) + -((V0O 2 + (V0 2 ) 2 ) + — (<pl + <j? 2 ) + (bi) 

|(0 4 + 4 ) + h{4>\ + 4) + £(§ 0102 + + 3^ 2 4 )) . 

For (3 = and 1, Eq.(Bl) corresponds to 0(l)xO(l) and 0(2) models, respectively. The 
evolution of b 2 for a "typical" value of G and H (we took G — 1, H — 3) is shown on Fig. 22. 
As it can be seen, the negative b 2 solution in 0(l)xO(l) never becomes positive as (3 gets 
increased and even more spurious solutions appear in the limit of 0(2). However, there 
is one-to-one correspondence between positive b 2 solutions in 0(l)xO(l) and 0(2). We 
found that the same is true for any other G and H . Thus, our most immediate exploration 
didn't realize the anticipated possibility that the negative b 2 solution in 0(1) x 0(1) theory 
corresponds to a positive b 2 solution in 0(2) theory. This may be a symptom of nonexistence 
of Goldstone mode in 1 + 1 dimension known as the Coleman's theorem [22]. Thus, further 
exploration including the higher dimension may be necessary for the better understanding 
of the spurious solution. 



4 In Ref. [9], Stevenson and Roditi used somewhat different notations for couplings. Their couplings 
a and (3 are related to our G and H as a = |p and (3 = t^t, respectively. The boundary of allowed 
a and (3 regions for the GEP were found as a curve (see Fig.l of Ref. [9]) that can be drawn by 
the following formulas: a = (2z + 3 + (z - 3)e z )/z 2 ,/3 = 3(z + 2 + (z - 2)e z )/z 3 with < z < oo. 
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FIG. 1. The nontrivial mass solutions for 
OR (solid and dashed lines), GEP (dotted 
line) and Hartree methods (squares corre- 
spond to the 4ttc 2 < 3 solution, circles to 
the 4-7TC 2 > 3 one). Here, G = g/m 2 and for 
G < 9.04 there are no nontrivial solutions. 
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FIG. 2. Vacuum energy densities of the 
nontrivial OR solutions (solid and dashed 
lines), GEP (dotted line) and Hartree meth- 
ods (squares correspond to the Attc 2 < 3 so- 
lution, circles to the Aire 2 > 3 one). Here, 
G = g/m 2 and for G < 9.04 there are no 
nontrivial solutions. 




-0.2 



FIG. 3. The effective classical potential corresponding to Fig. 1 with G = 8, 10 and 12. Al- 
though G = 8 has only a trivial effective potential denoted by a solid line, both G = 10 and 12 
have two nontrivial effective potentials denoted by I and II. 
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FIG. 4. Nontrivial solutions of the OR FIG. 5. Nontrivial solutions of the 

equations for (ft 4 and </> 6 theory, symmetric OR equations for 4 and <fi e theory, bro- 
phase. ken-symmetry phase. 
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FIG. 6. Nontrivial solutions of the OR equations for <p A and 6 theory, imaginary b phase. 



%,/m 2 




FIG. 7. Vacuum energy density of nontrivial solutions of the OR equations for 4 and (f> 6 theory, 
symmetric phase. The two plots (right and left) are identical but presented in two different angles 
of view. The 90° clockwise rotation of the left plot around the axis of £m/™ 2 is identical to the 
right plot. 
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FIG. 8. Vacuum energy density of non- 
trivial solutions of the OR equations for </> 4 
and </> 6 theory, broken-symmetry phase. 



FIG. 9. Vacuum energy density of non- 
trivial solutions of the OR equations for t/> 4 
and (f) 6 theory, imaginary b phase. 
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FIG. 10. Dimensionless couplings of the 
symmetric perturbative solution for G = 0. 
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FIG. 11. Dimensionless couplings of the 
broken-symmetry perturbative solution for 
G = 0. 
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FIG. 12. Dimensionless couplings of the 
broken-symmetry perturbative solution ver- 
sus the dimensionless coupling of the sym- 
metric solution (G = 0). 
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FIG. 13. Classical potential of bro- 
ken-symmetry (BS) and symmetric (S) the- 
ories at xf ] = 5-85 • 1(T 3 {H = 100, G = 0). 
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FIG. 14. Classical potential of bro- 
ken-symmetry (BS) and symmetric (S) the- 
1.33 -1(T 3 (H= 10 6 ,G = 0). 



, (6) 

ones at Xs 



FIG. 15. Classical potential of bro- 
ken-symmetry (BS) and symmetric (S) the- 
ories at xi 6) = 1.02- 10 -3 (H = 10 15 ,G = 0). 
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FIG. 17. The dimensionless effective potential V((f>) 
—(3 in notations of Ref. [9]), where 



V 



class 



+ £M 



/ m? for 



H/G = — 0.209(a = — f3 in notations of Ref. [9]), where <j) means ^m 88 • x axis is Re(4>), y 
axis is Im(4>) and z axis is V(<j>). In general, for the complex 4> the potential V(4>) is complex. For 
G =-0.838(solid line) and -1.68(long-dashed line), we show only the trivial solutions (t = and 
6 = 0) for the comparison with the nontrivial solutions for G = -2. 51 (short-dashed line),-3.35(dotted 
line) and -4.19(dot-dashed line). For each coupling constant, G = -2.51 (short-dashed line), 
-3.35(dotted line) and -4.19(dot-dashed line), there are two real b broken-symmetry solutions 
and one pure imaginary b solution. The two real b broken-symmetry solutions for each coupling 
constant(G =-2.51,-3.35 and -4.19) are shown on the Im((f>) = plane. For the pure imaginary 
b solutions, the potential graphs are drawn on the planes off from the Im(4>) = plane with the 
different imaginary b values depending on the G values. 
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FIG. 18. Real part of effective potentials 
for the spurious solutions of G=-0.838(solid 
line), -1.68 (long-dashed 

line), -2. 51 (short-dashed line), -3.35(dotted 
line) and -4.19(dot-dashed line). 
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FIG. 19. Imaginary part of effec- 
tive potentials for the spurious solutions 
of G=-0.838(solid line), -1.68 (long-dashed 
line), -2.51 (short-dashed line), -3.35(dotted 
line) and -4.19(dot-dashed line). 
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FIG. 20. Real part of effective potentials for the spurious solutions of G=-0.838(solid line), 
-1.68 (long-dashed line), -2.51 (short-dashed line), -3.35(dotted line) and -4.19(dot-dashed line) in 
the larger scale of classical fields. The blow-up of the results in a tiny region of 4>f^ ss — b < 0.025 
is shown in Fig. 18. 
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FIG. 21. The domain of trivial solutions 
in G, H parameter space (shadowed area). 
If the spurious solution shown in Fig. 6 is not 
taken into account, then the domain enlarges 
to the inside of the entire area bounded by 
the external solid line, which is equivalent to 
the GEP result shown in Fig.l of Ref. [9]. 
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FIG. 22. Correspondence between solu- 
tions of OR equations for 0(1) x 0(1) and 
0(2) (<p e )i +1 theories at G = 1.0, H = 3.0. 
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